A Mechanism for Polar Protein Localization in Bacteria 
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We investigate a mechanism for the polar localization of proteins in bacteria. We focus on the 
MinCD/DivIVA system regulating division site placement in the rod-shaped bacterium Bacillus 
subtihs. Our model relies on a combination of geometric effects and reaction-diffusion dynamics 
to direct proteins to both cell poles, where division is then blocked. We discuss similarities and 
differences with related division models in Escherichia coli and also develop extensions of the model 
to asymmetric polar protein localization. We propose that our mechanism for polar localization 
may be employed more widely in bacteria, especially in outgrowing spores, which do not possess 
any pre-existing polar division apparatus from prior division events. 



INTRODUCTION 

The polar localization of proteins is a vital component 
of the sophisticated internal organization of bacteria, and 
is important for many processes ranging from chemotaxis 
to accurate cell division El SUB- However, how such 
polar localization is achieved has remained largely un- 
known. The purpose of this report is to analyze one pos- 
sible mechanism for localization in a well-studied exper- 
imental context: the MinCD/DivIVA system regulating 
cell divisionpositioning in the Gram-positive bacterium 
B. subtilis |IllEllllaElElE3 We wiU then discuss 
how these ideas might be applied more generally to the 
generation of polar localization/asymmetry in bacteria. 

The dynamics of bacterial cell division have been in- 
tensively studied, particu larly in B. subtilis la la B B, 

iHTnUp and e. coli El El El El El El Ea lia 

2 It I22L |2jj, where the polar localization of proteins is 
known to be crucial. In both cases, the location of the 
division site is determined by the tubulin-like GTPase 
FtsZ which assembles into a ring-shaped structure just 
under the cell membrane at midcell. In E. coli accurate 
targeting of the FtsZ ring to midcell is achieved both 
by the inhibitory effect of the nucleoids ("nucleoid oc- 
clusion") and through the MinCDE proteins MinC 
and MinD form a complex on the cell membrane, which 
undergoes coherent pole-to-pole oscillations, with a pe- 
riod of around a minute, driven by MinE. As a result 
of the oscillations, the time-averaged MinCD concentra- 
tion is minimized at midcell and maximized at the cell 
poles. Since MinC inhibits FtsZ polymerization, the FtsZ 
ring assembles at midcell where the average MinC con- 
centration is lowest. In this way, through the temporary 
sequestration of MinCD at the cell poles, the midcell can 
be accurately targeted for cell division. 

Recently, several groups have employed simple 
reaction-diffusion mathematical models (either deter- 
ministic HI HI or stochastic H3) to try to un- 
derstand the MinCDE dynamics (see also Refs. [H 
for related models). The principal result of this model- 
ing effort is that the oscillations are a spontaneous self- 



organized phenomena arising from an intrinsic instabil- 
ity of the reaction-diffusion dynamics. A vital element 
of this instability is the varying diffusion constants of 
the MinCDE proteins: high in the cytoplasm; low when 
polymerically bound to the membrane. Crucially, the 
oscillations form spontaneously as a result of this dispar- 
ity of the diffusion constants and of the nonlinear bind- 
ing/unbinding interactions between the proteins. As a 
result, the oscillatory dynamics does not require any pre- 
existing topological markers to distinguish various loca- 
tions within the cell, as had previously been supposed. 

Very recently, improved fluorescent imaging has re- 
vealed that the E. coli Min proteins polymerize into sta- 
ble helices on the cytoplasmic membrane |22j . This phe- 
nomenon is not explicitly included in any of the existing 
mathematical models, which may therefore require some 
amendment. Nevertheless, the models generally agree 
well with fluorescent imaging experiments, especially in 
the observation of periodic MinCDE/FtsZ bands in fil- 
amentous mutants. These results are very difficult to 
understand other than as a pattern-forming reaction- 
diffusion instability. Hence, we believe the essentials of a 
reaction-diffusion mechanism underpinning the MinCDE 
oscillations are likely to be correct. Although the forma- 
tion of the Min proteins into filaments has only been ob- 
served in E. coli (and also in vitro H^), this observation 
may have important implications in other bacteria. Al- 
ready in B. subtilis the shape-determining proteins MreB 
and Mbl are known to form helical cables that could com- 
prise cytoskeletonal-like elements 

Despite possessing two of the three same proteins 
(MinC and MinD, but not MinE), B. subtilis appears to 
use a quite different mechanism to achieve central place- 
ment of the FtsZ ring in vegetative cells H H H H 
E1EJE3- In B - subtilis, MinCD do not oscillate, but 
instead are localized to both cell poles, where they are 
anchored by the protein DivIVA (which is structurally 
unrelated to MinE). Also involved are the (initially) ho- 
mogeneously distributed proteins ZapA, which promotes 
|3l| . and EzrA, which inhibits [22]> FtsZ ring assembly. 
Quite how DivIVA is able to recognize the cell poles 
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has remained unresolved and forms the main object of 
this report. With MinCD localized to the cell poles, the 
FtsZ ring is again unable to nucleate there. However, 
it is likely that additional mechanisms are used to ac- 
curately locate the FtsZ ring to midcell |33|, since the 
MinCD/DivIVA system functions primarily to block po- 
lar division These additional (and so far unchar- 
acterized) mechanisms may be especially important as 
nucleoid occlusion is ineffective as a means of positioning 
in B. subtilis [Sfij . Despite this, the precision of midcell 
division is the same in E. coli and B. subtilis |34j |. even 
though their regulatory mechanisms differ considerably. 

Polar protein localization is not only restricted to po- 
sitioning of the division site. For example in Caulobac- 
ter crescentus, polar asymmetry is exploited to ensure 
that daughter cells are dissimilar: one a swarmer cell 
with a polar flagellum, and one a stalked cell with a po- 
lar stalk that can adhere to a source of food 0] . In 
each case, the development of a cell with a polar stalk 
or flagellum is critically dependent on spatiotemporally 
controlled protein localization, and on the subsequent 
proteolysis and release of these components at the ap- 
propriate time. Other examples of polar localization in- 
clude thepolar clustering of chemoreceptors, in both E. 
coli H^,|33, Hl| and C. crescentus [3^, and also the po- 
lar localization of the virulence protein IcsA ^3], which 
promotes undirectional actin-based motility of the bac- 
terium Shigella flexneri. 



POLAR PROTEIN LOCALIZATION IN B. 
SUBTILIS 

Polar division in vegetative cells of B. subtilis is blocked 
by the presence of membrane localized MinC. MinC is 
normally a cytoplasmically localized protein, however, 
in the presence of MinD the two interact strongly and 
bind to the membrane ,8|. Polar localization is facili- 
tated by the action of the protein DivIVA, which acts to 
anchor MinCD to the poles 0Q. Targeting of MinCD 
to the poles, is dependent on DivIVA, as in filamentous 
cells lacking DivIVA, MinCD was partly distributed dif- 
fusively along the cell 00- Furthermore the localization 
and activity of MinD is believed to require ATP binding 
and possibly hydrolysis Q. However, MinCD can also 
(MinC more strongly than MinD) target past/present 
division locations independently of DivIVA E| ■ Simi- 
larly in a minD mutant, DivIVA simply localizes to divi- 
sion sites regardless of location, again demonstrating an 
affinity for some component of the cell division appara- 
tus 0. Thus in wild-type cells, MinCD/DivIVA local- 
ize to the poles (old division sites), blocking cell division 
there. However, before septal constriction occurs and the 
cell divides, MinCD/DivIVA assemble approximately si- 
multaneously at midcell 0] , and are retained there after 
division [g. In this way MinCD is delivered to the 



new poles of the two daughter cells, meaning that polar 
division is immediately prevented at the new poles. 

This observation immediately provides one possible ex- 
planation for polar localization: since the cell division 
apparatus marks the location of the new cell poles in 
the daughter cells, if the relevant proteins can localize 
to a component of the division apparatus, then the cell 
pole will continue to be marked in the daughter cells. 
However, one may question whether this mechanism is 
really sufficiently robust: if MinCD/DivIVA fail to cor- 
rectly localize to the site of division, then the poles of the 
resulting daughter cells will remain unmarked, meaning 
that division in those cells would be corrupted. It would 
certainly therefore be desirable to have other means of 
identifying the cell poles. Experiments in outgrowing 
spores of B. subtilis indicate that this is indeed the case 
[Hi D3 MinCD/DivIVA are targeted to the cell poles 
very rapidly after spore germination, even in cells lacking 
FtsZ. Thus, polar localization can occur even in germi- 
nating spores which do not possess pre-existing division 
apparatus left over from prior division events. 

So the question arises, what other mechanism could 
be responsible for the polar localization? One possible 
answer is that a self-organized reaction-diffusion insta- 
bility exists, similar to that already studied in E. coli. 
Of course, for B. subtilis, the instability would have to be 
to a steady-state distribution where the MinCD/DivIVA 
do not oscillate (i.e. in linear stability analysis one would 
expect a real rather than a complex growth eigenvalue). 
With appropriate interactions and diffusion constants, 
such an instability could certainly exist. However, one 
of the key predictions of such an instability is the exis- 
tence of a characteristic wavelength, depending on the 
microscopic parameters of the model, but independent 
of the geometry 0, HE E§] ■ Thus in long filamentous 
cells of E. coli, several wavelengths can be fit into a 
single cell and one therefore expects to see multiple os- 
cillating bands. The observation of this phenomena in 
filamentous E. coli 13] is powerful evidence in favor of 
the reaction-diffusion mechanism. The situation in B. 
subtilis is, however, rather different: in filamentous cells 
lacking FtsZ, MinCD/DivIVA are primarily localized to 
the cell poles 0] , without any regularly spaced inter- 
mediate bands in between, as one would expect for a 
reaction-diffusion instability. One could argue that the 
affinity of MinCD/DivIVA for the pre-existing (polar) 
division apparatus from previous division events might 
mask the existence of any intermediate bands. However, 
in outgrowing spores, which have no pre-existing division 
apparatus, there are still no intermediate bands, even 
in long filamentous mutants (up to about 13 fim long) 
[ill Ha ], Thus we see no evidence for the existence of 
a characteristic wavelength. Of course, it is just possi- 
ble that the characteristic wavelength is simply very long, 
too long to have been seen in past experiments; but, there 
is currently no hint of this. 



3 



Instead we propose that it is a geometric effect in com- 
bination with reaction-diffusion dynamics which marks 
the end of the cell. The essence of our approach is to as- 
sume that MinCD has a reduced rate of membrane bind- 
ing at the cell poles. If, for example, MinCD polymerize 
into filaments, as in E. coli 23], this might be due to 
the increased curvature of the cell poles interfering with 
membrane polymerization. Since DivIVA may then bind 
to the edges of MinCD clusters and so stabilize them, 
the MinCD (and DivIVA) concentrations end up being 
maximized close to the cell poles, even though MinCD 
initially has a reduced rate of binding there. Similarly, 
when the cell divides a region of high curvature is created 
at midcell, again attracting the MinCD /DivIVA proteins 
independent of the division apparatus. In the next sec- 
tion, we define this model more precisely and report on 
the results of numerical simulations. 



THE MODEL 

Previous modeling work on the MinCDE system has 
demonstrated the effectiveness of simple deterministic 
reaction-diffusion equations |2J, |2jJ |2y| . Accordingly we 
follow a similar modeling approach here. The dynamics 
of MinC is similar to that of MinD, so, for simplicity, we 
drop MinC from our model. The fact that MinD can lo- 
calize normally in minC mutant cells also supports this 
approach ||. Initially, we also assume a "conservative" 
model, where the MinD/DivIVA are able to bind/unbind 
from the membrane but are not degraded in the process. 
This is again similar to some of the models used for E. 
coli 0> O E3- The relevant set of one-dimensional 
reaction-diffusion equations are: 
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Here, the first equation is for the cytoplasmic MinD den- 
sity pD, the second for the membrane bound MinD den- 
sity pd, the third for the cytoplasmic DivIVA density pa, 
and the fourth for the membrane bound DivIVA density 
p a . In the above model, we assume that MinD binds to 
the membrane with rate Ad, in agreement with green flu- 
orescent protein (GFP) fusions which indicate that MinD 
is a membrane-localized protein 0, 0- However, one 
would also expect that the rate of binding is reduced as 
the MinD membrane density increases. This accounts 



for the X' D term. We also assume that MinD can unbind 
spontaneously (Ad term), but that this is suppressed by 
DivIVA (X' d term) 0, @- It is also possible to insert a 
factor (1 + A^'pd) -1 in this term, which would model in- 
creased stabilization of membrane MinD from clustering. 
However, we have found that our results were not sig- 
nificantly affected by this term, which we have therefore 
left out for simplicity. Notice that the DivIVA functions 
to stabilize the membrane MinCD, hence its role is very 
different from the E. coli MinE protein which destabilizes 
membrane MinCD. This difference is partly responsible 
for the structure of the equations being rather different 
in this case from the earlier models for E. coli MinCDE 

HI IH ■ DivIVA is also believed to be an intrinsi- 
cally cytoplasmic protein [jj. Hence, in order to bind to 
the membrane in outgrowing spores (without pre- 
existing division apparatus), it presumably requires the 
presence of MinD. We assume here, similar to Ref. HJ, 
that it binds to (and then stabilizes) the edges of MinD 
clusters. This could correspond to DivIVA "capping" the 
ends of MinCD filaments 91. This process would then ac- 
count for the Xa gradient term \dpd/dx\, which is max- 
imized at the edges of a cluster. Note that the gradient 
term has a modulus sign, as the rate of binding should 
not depend on the sign of the gradient. Again we assume 
that the DivIVA binding is suppressed as the membrane 
DivIVA density increases (A^ term). Finally DivIVA can 
spontaneously unbind (A a term), but this is suppressed 
by DivIVA itself (X' a term). This indicates that DivIVA 
itself oligomerizes/forms clusters for maximum stability, 
for which there is also experimental evidence 

The above assumptions are fairly minimal and are all 
closely based on the available experimental evidence of 
GFP fusions and deletion mutants 0, 0, S H ES [H 
H^ |. As we have said, we then add one extra ingredient, 
namely a position dependent membrane binding Xr>(x) 
for the MinD, where Xd(x) is reduced in regions close to 
the cell ends (see below). 

For our computer simulations, we assume hard wall 
boundary conditions (i.e. dpi/dx = 0, i — {D,d, A, a}, 
at the boundaries) and we numerically integrate the 
above equations, with lattice size dx — 0.02 pm and time 
increment dt = 1 x 10 -5 s. We start from random initial 
conditions, where the MinD and DivIVA are randomly 
distributed throughout the cytoplasm. We begin the sim- 
ulation with a total of 50 lattice sites (bacterial length 

1 pm), with initial densities of 1000 pmT 1 for MinD and 
750 pmT 1 for DivIVA. The actual concentration levels 
for B. subtilis have not been accurately measured, so we 
assume concentration levels similar to those observed in 
E. coli |4l|. During the course of the simulation, we 
then allow the bacterial length to increase, modeling a 
growing bacterium, with growth rate of 4 pm /hour. In 
other words, every 1.8 x 10 5 time steps an extra site is 
inserted, randomly located. Note, however, that the ex- 
tra site is never inserted at the inert cell ends 01 (de- 
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fined here by regions where Xo{%) is reduced, see below). 
As the bacteria grows in the simulation we add extra 
MinD/DivIVA, randomly localized in the cytoplasm, to 
ensure that the average densities remain constant (i.e. 
the total number of MinD/DivIVA proteins is propor- 
tional to the cell length). There is currently no informa- 
tion on whether the average densities do actually remain 
constant as the cell grows, but the assumption used here 
is the simplest possible. Nevertheless, precise regulation 
of the densities turns out not to be particularly impor- 
tant. Considerable overproduction eventually leads to 
the MinCD/DivIVA proteins saturating the poles and 
spilling over into the reminder of the cell length. Un- 
derproduction simply translates to less MinCD/DivIVA 
being present along the cell length in between the poles. 
The important point is that the bacterium has to manu- 
facture sufficient MinCD/DivIVA to ensure that the cell 
poles are always well populated with these proteins, but 
not so much that the proteins spill over to midcell lead- 
ing to division inhibition. Note that in our simulations, 
we exclusively model the case of outgrowing spores with- 
out pre-existing division apparatus. Clearly, modeling 
protein localization to pre-existing binding sites in wild- 
type vegetative B. subtilis cells is straightforward, and 
not particularly revealing. 



In our simulations we use the following parameter 
set: D D = 2.0 /xr^s" 1 , D d = 0, D A = 2.0 /im 2 s -1 , 
D a = 2 x 10~ 4 imps- 1 , X D = 10 s" 1 , X' D = 0.001 fim, 
X d = 3.0 s" 1 , X' d = 0.2 fim, X A = 0.03 /WV 1 , 
A^=8x 10~ 4 fim, X a = 2.0 s"\ X' a = 0.005 fim. Fur- 
thermore, for the closest 10 sites to both ends of the 
bacteria, the binding rate A_d is reduced by a factor of 
2 (our results are robust to changes in this suppression 
factor). It is also possible to generalize the model so 
that the MinD /DivIVA proteins are degraded at a (slow) 
rate, whilst increasing the production rate to keep the 
average densities constant. We have verified that this 
modification also makes no qualitative difference to our 
results. The values for the above cytoplasmic diffusion 
constants are of the same order as that measured (for 
different proteins) in the cytoplasm of E. coli |43J . As for 
the membrane diffusion constants, we have used similarly 
small values to those used in E. coli [3 E3 • Previously 
this was an important component in the spontaneous 
reaction-diffusion instability; this is, however, not the 
case here: the small membrane diffusivities are now nec- 
essary to preserve the polar density concentration max- 
ima and to prevent the concentrations from equalizing 
everywhere along the cell membrane. For the remaining 
parameters there is no experimental guidance, but we 
have taken care to check that our results are robust to 
variations in these parameters. 



RESULTS: POLAR PROTEIN LOCALIZATION 

In Figure 1, we show the result of our computer sim- 
ulation of the spore outgrowth experiment. We see 
very similar results to those observed experimentally 
Starting from random initial conditions, the 
MinD/DivIVA rapidly (within a few minutes) localize to 
the poles of the cell. This agrees with experiments indi- 
cating that MinD/DivIVA localize to the cell poles very 
early in spore germination [ill Il2| . In our model the 
mechanism for this localization follows from the reduced 
rate of MinD binding to the cell ends. This reduced polar 
value for Xo(x) is absolutely critical for polar localiza- 
tion as it sets up a pronounced membrane MinD density 
gradient close to the cell poles. The DivIVA then binds 
to these positions of increased gradient, stabilizing the 
MinD/DivIVA. In this way the MinD/DivIVA concen- 
trations end up being maximized close to the cell ends. 
In fact, the concentrations are not maximized precisely 
at the poles, where there is actually a dip (see Figures 2 
and 3). However the polar concentrations still end up 
being much higher than the densities well away from the 
cell ends. As the cell continues to grow, the polar regions 
become saturated with MinD/DivIVA and the newly pro- 
duced proteins cannot localize there. Instead the pro- 
teins localize to random membrane patches along the cell 
length. Furthermore, these patches appear to propagate 
along the cell length (see Figure 1). This movement re- 
sults from the coupling of the DivIVA to the edges of 
MinD clusters, slightly away from the local MinD den- 
sity maxima. Since MinD is stabilized by the DivIVA this 
encourages MinD to rebind slightly away from its current 
density maxima. This in turn shifts the position of the 
maximum MinD density gradient, thus also encouraging 
the DivIVA to rebind in a shifted location. In this way 
the coupled MinD /DivIVA density maxima both tend to 
propagate along the membrane, as shown in Figure 1. It 
would be interesting to see if such movement can be ob- 
served in experiments. A snapshot of these patterns can 
be seen in Figure 2, where we show the MinD/DivIVA 
membrane density profiles at length 10 fim. The DivIVA 
pattern is very similar to that observed experimentally in 
filamentous spore outgrowth experiments 0, , with 
sharp polar peaks and less intense randomly distributed 
membrane patches along the cell length. Note that this 
disordered pattern is seen everywhere inside the simu- 
lated outgrowing spore (including close to the cell ends) 
if Xo(x) is a constant, emphasizing again that polar pro- 
tein localization is critically dependent on a position- 
dependent MinD binding rate. 

More quantitative data, with profiles of the membrane 
MinD/DivIVA densities at different times are shown in 
Figure 3(a,c,e). In Figure 3(b,d,f) we show examples of 
the membrane density averaged over 90 s as a function of 
position for various total cell lengths. A 90 s average is 
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FIG. 1: Spacetime plot of membrane-bound DivIVA (left) 
and MinD (right) densities in a simulated outgrowing spore. 
Time increases from top to bottom and the simulation lasts 
8100 s, starting with length 1 fim and ending with length 
10 fim. White (black) corresponds to areas of high (low) 
concentration. 




FIG. 2: Snapshot of the DivIVA (top) and MinD (bottom) 
membrane density profiles at length 10 fim. White (black) 
corresponds to areas of high (low) concentration. Note that 
these are one dimensional simulation results (two dimensional 
presentation in the figure is for clarity). 

chosen as, at least in E. coli, this is the rough time over 
which the FtsZ averages in locating the FtsZ ring [lj| . As 
can be seen, the MinD concentration is nicely minimized 
at the cell center for shorter cell lengths. However, as the 
cell length increases, the precision with which the center 
can be identified is considerably reduced, with the aver- 
aged signal fluctuating considerably at 8 /jm. This data 
is consistent with evidence that the MinCD /DivIVA sys- 
tem is primarily responsible for blocking polar division, 
rather than in accurately placing midcell division |34| . 
Notice also how the MinD is less tightly localized at the 
cell poles; its signal decreases more slowly than for the 
DivIVA towards the center of the cell. This also agrees 
with experiment Q. 

RESULTS: ASYMMETRIC PROTEIN 
LOCALIZATION 

Up to this point, we have modeled symmetric protein 
localization, such as that seen in the MinCD/DivIVA sys- 
tem in B. subtilis. However, it is also interesting to ask 
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FIG. 3: (a,c,e) Simulations of outgrowing spores: instanta- 
neous membrane MinD (dashed line) , and membrane DivIVA 
(full line) density profiles at lengths 2 fim, 4 /jm and 8fim; 
(b,d,f) Membrane MinD (dashed line) and membrane DivIVA 
(full line) density profiles averaged over 90 s at lengths 2 fim, 
4 fim and 8 fim. Maximum densities for each of the membrane 
MinD and DivIVA normalized to unity. 

whether one can modify the model to obtain asymmet- 
ric localization, i.e. to one pole only. This modifica- 
tion turns out to be quite straightforward: in essence 
one merely has to remove the X' D and A^ terms which 
restrict membrane binding when the membrane densi- 
ties are already high. In this case, starting from ran- 
dom initial conditions, the dynamics rapidly localizes 
the proteins to both cell poles. However, one randomly 
chosen pole will, due to fluctuations in the random ini- 
tial conditions, gain a slightly higher density than the 
other. Since higher densities yield greater stability for 
the membrane aggregations (due to the A^ and X' a terms) , 
and since there are no longer saturation effects, the fa- 
vored pole, with an initially slightly higher concentra- 
tion, will continue to attract additional proteins from 
the other pole, via cytoplasmic diffusion and then po- 
lar membrane capture (cf. Ref. 01)- Of course this 
process will also occur in reverse from the favored to the 
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less favored pole, but because of its higher density, the 
net rate of accumulation will be towards the pole with 
the higher initial concentration. In this way, the less fa- 
vored pole is eventually drained of its protein population. 
Starting from randomly distributed cytoplasmic MinD 
and DivIVA, with average densities 1000 [imT 1 , and 
750 [imT 1 , respectively, an appropriate parameter set is: 
D D = 2.0 /iTO 2 s-\ D d = 0, D A = 2.0 /tmV 1 , D a = 
2 x 10~ 4 ^m 2 s _1 , \ D = 10 a -1 , X' D = 0, X d = 3.0 s _1 , 
A^ = 0.001 ^m, A A = 0.012 ^m 2 s _1 , A^ = 0, A a = 
1.0 s _1 , AJj = 1.2 x 10~ 4 fim. In this case localization to 
a (randomly selected) pole occurs in about 800 s (data 
not shown). Alternatively, replacing some of the above 
parameters with A^ = 1 x 10 -4 /jm, Aa = 0.002 ^m 2 s _1 , 
A a = 1.6 s , and X' a = 1.6 x 10~ 4 /xm, unipolar localiza- 
tion occurs far more rapidly, within about 40 s. A process 
analogous to this may be happening in asymmetric divi- 
sion in sporulating B. subtilis 0] (see below), although 
in that case the protein SpoIIE may also play a crucial 
role 0. 

DISCUSSION 

In this paper, we have proposed a mechanism for the 
symmetric or asymmetric polar localization of proteins, 
relying on a combination of geometric information and 
reaction-diffusion dynamics. This mechanism could be 
particularly important for polar protein localization in 
cells growing out from spores, where there are no pre- 
existing division apparatus from prior division events. 
However, even in vegetative cells, where the poles could 
be marked by the division apparatus from prior divisions, 
the proposed mechanism could be useful for increased re- 
liability (and redundancy) in locating the poles. We have 
focused particular attention on the MinCD/DivIVA sys- 
tem in B. subtilis, where our results agree well with exper- 
iments on outgrowing spores. We have also made a num- 
ber of testable assumptions: very low diffusion constants 
for membrane-bound MinCD/DivIVA; DivIVA binds to 
the edges of MinD clusters; and membrane binding of 
MinD is heavily influenced by geometric effects (possi- 
bly by the high curvature of the cell poles restricting 
MinD membrane polymerization). We should point out, 
however, that experiments in round cell mutants have 
indicated that simple curvature of the cell membrane is 
insufficient to elicit DivIVA targeting [47^ . Nevertheless, 
this does not contradict our results: according to our 
model, what is important is that there be a difference 
in curvature: one part of the cell with low curvature, 
another part (the cell poles) with high curvature. 

As we mentioned in the introduction there are several 
other regulatory systems in bacteria with similar polar 
localization properties. For exam ple, the MCP chemo- 
taxis receptor proteins (in E. coli) |'id I37I l3S| and McpA 
(in C. crescentus) [39| are both polar localized. Fur- 



thermore, filamentous mutants in both cases show an 
intriguing pattern of localization: the proteins are still 
polar localized, but in addition there is a seemingly ran- 
domly distributed component along the filament length 
between the two poles. This pattern is qualitatively rem- 
iniscent of the pattern seen (and modeled in this paper) 
for MinCD/DivIVA in filamentous B. subtilis mutants 
grown from spores 0, 0] . Although, for the chemotaxis 
cases, it may simply be that the proteins have saturated 
polar binding sites and arc then "spilling" over into the 
remainder of the cell. However it is again possible that 
there may be a geometric component to the localization, 
similar to that proposed in this report. This may be even 
more likely in the case of chemotactic proteins in out- 
growing spores. Already the chemotaxis protein McpB is 
known to be polar localized in vegetative B. subtilis ^3 • 
It would be very interesting to study its localization prop- 
erties in outgrowing spores without pre-existing division 
apparatus. In addition, it has been shown that DivIVA 
localizes to division sites both in E. coli and, unexpect- 
edly, in the fission yeast Schizosaccharomyces pombe ^3] • 
However, this may again be due to binding to division 
apparatus proteins (possibly actin in S. pombe). 

Given that the geometric properties (high curvature) 
of the cell poles may be important in influencing MinD 
binding, it is, of course, important to understand how 
the cell shape is maintained. In B. subtilis, cell elonga- 
tion and maintenance of the cell cylinder is achieved by 
Mbl-directed dispersed helical insertion of cell wall mate- 
rial throughout the cylinder, while the hemispherical cell 
poles are effectively inert |42j. In vegetative cells, where 
a midcell septum is formed, it is plausible to see how 
the hemispherical end is constructed and thus how over- 
all cell shape maintained. However, in spore outgrowth, 
this is not so obvious. B. subtilis spores themselves are 
already ellipsoidally shaped, so the initial state for spore 
germination already has the ends marked via high cur- 
vature. However quite how this shape is produced in the 
earlier sporulation event is rather unclear. 

It is also intriguing to note that, although MinC and 
MinD are used by both B. subtilis and E. coli, the two 
bacteria appear to use entirely different physical mecha- 
nisms to exclude division from polar sites. E. coli uses a 
spontaneous reaction-diffusion instability, whereas geo- 
metric effects (reduced polar binding of MinD) and asso- 
ciated reaction-diffusion kinetics, appear to be employed 
in B. subtilis. Interestingly, it is actually a reduced rate 
of MinD membrane binding at the cell poles which, in 
our model, eventually leads to polar localization due to 
DivIVA binding/stabilizing the edges of MinD membrane 
clusters. It would be interesting to further study these 
mechanisms with their evolutionary origin in mind, es- 
pecially since both employ reaction-diffusion dynamics. 
Although we believe that the mechanism outlined in this 
paper provides a very plausible scenario for polar protein 
localization, other mechanisms are still possible. For ex- 
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ample, one could imagine that the DivIVA makes the 
poles more attractive for MinCD polymerization, per- 
haps by altering the structure of the MinCD filaments 
so that they prefer a higher curvature region (see also 
Ref. 01 )■ Given that the model described in this paper 
provides several testable predictions, experiments should 
be readily able to discriminate between these possible 
models. 

One may also speculate on why E. coli and B. subtilis 
have differing mechanisms. One possible explanation is 
that B. subtilis must have extra flexibility in its cell di- 
vision regulation as, unlike E. coli, it sporulates when 
confronted with harsh environmental conditions. In the 
process of forming a spore, the division septum is relo- 
cated from midcell to one of the cell poles. Hence B. 
subtilis must have the ability to undo its polar division 
block. All that would be required is for another protein 
to bind to DivIVA, displacing the MinCD, thus freeing 
up the polar sites, as suggested in Ref. |49|. A suitable 
candidate for such a displacing protein has recently been 
shown to be RacA j^. The MinCD presumably then 
relocates to the cell center, forcing the FtsZ ring to spi- 
ral outwards to the cell poles, as observed in Ref. |45| . 
As only a single pole is then selected for spore formation, 
the cytokinetic proteins should localize asymmetrically, 
possibly in a similar way to that discussed above. 

Finally, we mention that the repertoire of division site 
positioning is likely still not exhausted by the MinCDE 
and DivIV A sy stems analyzed here and in previous mod- 
eling work 123, |26j, |2jJ . C. crescentus is able to divide 
asymmetrically with high precision even though it lacks 
all of the MinCDE and DivIVA proteins . How this is 
achieved remains unknown and hints that there may be 
yet other mechanisms capable of providing accurate po- 
sitional information and hence of regulating subcellular 
architecture. 
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